Quantum synthesis of arbitrary unitary operators 
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Nature provides us with a restricted set of microscopic interactions. The question is whether we 
can synthesize out of these fundamental interactions an arbitrary unitary operator. In this paper we 
present a constructive algorithm for realization of any unitary operator which acts on a (truncated) 
Hilbert space of a single bosonic mode. In particular, we consider a physical implementation of 
unitary transformations acting on 1-dimensional vibrational states of a trapped ion. As an example 
we present an algorithm which realizes the discrete Fourier transform. 
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I. INTRODUCTION 

Controlled manipulations with individual quantum 
systems such as trapped ions or cold atoms in atomic 
physics, molecules irradiated by laser fields, and Ryd- 
berg atoms interacting with quantized micromaser fields, 
provides us with a deeper understanding of fundamen- 
tal principles of physics. Simultaneously, the possibility 
to control individual quantum systems opens new per- 
spectives in application of quantum physics. Specifically, 
coherent control over dynamics of quantum systems is of 
vital importance for quantum computing and informa- 
tion processing |l],|2| . 

Quantum information processing can be schcmatici- 
cally devided into three stages. The first stage is the 
encoding of information into quantum systems, i.e. this 
corresponds to a preparation of states of quantum sys- 
tems. The second stage is the information processing 
which in general is equivalent to a specific unitary evo- 
lution of the quantum system, i.e. this is an application 
of a given quantum algorithm. The third stage is the 
reading of output states of quantum registers (i.e. the 
"decoding" of information from quantum systems). Ob- 
viously this final stage is the measurement of quantum 
system and the reconstruction of relevant information. 

There are several physical systems which are believed 
to be hot candidates for quantum processors. In partic- 
ular, Cirac and Zoller have shown that a system of 
trapped ions can be utilized as a prototype of a quantum 
computer. Therefore it is of great interest to understand 
how the three stages of the information processing as 
specified above can be implemented in this system. 

(i) State preparation: Recently, several methods for de- 
terministic synthesis (preparation) of vibrational states 
of trapped ions have been proposed. In particular, a 
scheme for preparation of quantum states of one and two- 
mode bosonic fields (e.g. 1 and 2-D quantum states of 
vibrational motion of trapped ions) have been proposed 
by Law and Eberly Q, and Kneer and Law J5j (see also 
H and for a more general discussion on the state prepa- 
ration see 0). 

(ii) State measurement: There are various experimen- 



tal techniques which allow to measure and reconstruct 
quantum states of trapped ions (for a review see S). 

(Hi) Arbitrary unitary evolution: One of the most im- 
portant task in information processing is to design pro- 
cessors which take an arbitrary input and process it ac- 
cording to a specific prescription. Nature provides us 
with a restricted set of microscopic interactions. The 
question is whether we can synthesize out of these fun- 
damental interactions an arbitrary unitary operator. In 
this paper we present a constructive algorithm for real- 
ization of any unitary operator which acts on a finite- 
dimensional Hilbert space. An algorithmic proof that 
any discrete finite-dimensional unitary matrix can be fac- 
torized into a sequence of two-dimensional beam splitter 
transformations was given by Reck et al. |J. The prob- 
lem of controlled dynamics of quantum systems has been 
addressed recently by Harel and Akulin [JlO| and by Lloyd 
and Braunstein Harel and Akulin |(l0|| have pro- 

posed a method to attain any desired unitary evolution 
of quantum systems by switching on and off alternatively 
two distinct constant perturbations. The power of the 
method was shown in controlling the 1-D translational 
motion of a cold atom. 

Our aim is to find a constructive algorithm to realize 
an arbitrary unitary operator V which transforms any 
state \ip) of a single bosonic mode, e.g., an 1-D vibrational 
state of a trapped ion in x direction, to another state \yj'}, 
i.e., \ip'} = V\ip). In particular, we consider a truncated 
(N + l)-dimensional Hilbert space of the bosonic mode. 
Within this truncated Hilbert space the desired unitary 
operator V in the number-state basis reads 
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Under action of the operator V a given input state trans- 
forms as 
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where c' n = X)m=o Vn,mC-m- Our task is to represent by a 
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feasible physical process any operator V (|I|) . The synthe- 
sis of operators enables thus to realize universal quantum 
gates for qubits which can be encoded into vibrational 
levels. 

The paper is organized as follows. In Section |J we 
briefly introduce physical tools which we use to realize 
arbitrary unitary operators for a trapp ed i on. The syn- 
thesis algorithm is described in Section III . The method 
is illustrated in Section IV where a realization of quantum 



gates which perform the Fourier transform is considered. 
We also discuss stability of the algorithm and possibili- 
ties to realize also non- unitary operators. We finish our 
paper with conclusions. 



II. TOOLS FOR SYNTHESIS: LASER 
STIMULATED PROCESSES 

Our realization of the unitary transformation V given 
by Eq.(Q) is based on an enlargement of the Hilbert sub- 
space of the given system. Namely, the transformed 
bosonic mode corresponds to one vibrational mode (in 
x direction, for concreteness) of a quantized center-of- 
mass motion of an ion confined in the 2-D trapping po- 
tential. Within our synthesis procedure the vibrational 
x-mode becomes entangled with the auxiliary degrees of 
freedom (ancilla) which are represented by the second vi- 
brational mode (e.g., quantized vibrational motion in y 
direction) and three internal electronic levels \a), \b), \c) 
of the ion. The particular choice of the physical system is 
motivated by a feasibility of highly coherent control over 
motional degrees of freedom as demonstrated in recent 
experiments || which makes trapped ions to be a hot 
candidates for quantum processors. 

A physical realization of the desired operator V for an 
ion confined in a 2-D trapping potential consists of a se- 
quential switching (on/off) of laser fields which irradiate 
the ion. Namely, we utilize four types of laser stimulated 
interactions which are associated with the following (ef- 
fective) interaction Hamiltonians: 



H^ m ) = {A y +s x ){\a){a\ 
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where 
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and T{a q a q ) = e V J2 k =o ( k +iy.k\ 

Dicke parameters are denned as r\ q = oj q / {c^2m a v q ) 
(assuming units such that h — 1) where uj q ,v q are fre- 
quencies of the laser and vibrational mode in direction q 
(q = x, y), respectively. Further, lP m denote the Laguerre 
polynomial; m a is the mass of the ion. 

The dynamical Stark shift operator 
Sx = J2 m s m\m) x (m\ is induced by a detuned standing- 
wave laser field applied in x direction ||] . For x 9i the 



effective Hamiltonian H^ 1,m ' addresses the states with 
fixed number m of excitations (indicated by the super- 
script) in the mode x setting the detuning A y of the laser 
field, applied in the orthogonal y direction, equal to the 
dynamical Stark shift — s m . In other words, when the 
interaction is governed by the Hamiltonian H^ 1,m ' with 
A y = — s m there is an exchange of the population only 
between the states |m, n) ® \a) <^> \m, n) <£> \b) with the 
given number of quanta m in the vibrational mode x and 
any number of phonons n in the mode y. The populations 
of the other number states (with number of excitations 
in the mode x different from to) effectively do not change 
due to a large detuning. However, there are significant 
phase shifts of the amplitudes of the off-resonant states. 
It should be stressed that the addressing of states with a 
given number to of excitations in the mode x is effective 
only for large ratios x/<?i ^ !• Moreover, the approx- 
imate Hamiltonian ifC 1 '™) is itself justified only when 
X ^ I Si I (f° r details see 

The considered interactions ([}]) are quite typical for a 
trapped ion. The effective interaction Hamiltonians rep- 
resent a classical driving of the ion [see H^] and a non- 
linear Jaynes-Cummings model [12| [see iJ( 2 )>( 4 )] when 
the lasers applied in given directions (x,y) are tuned to 
appropriate vibrational sidebands. These Hamiltonians 
are thoroughly discussed in papers |^,^[ . 

The dynamics governed by the interaction Hamiltoni- 
ans (H) can be separated into independent 2-D subspaces. 
Switching on a particular interaction "channel" associ- 
ated with one of the interaction Hamiltonians H^] (||) 
for a time r is described as the action of the correspond- 
ing unitary time-evolution operator U^> = exp[— 
on the state vector of the system under consideration. 



III. SYNTHESIS OF TRANSFORMATIONS 

To implement the desired transformation V (Q) for an 
ion confined in the 2-D trapping potential we realize a 

mapping \ip x , y ) — >|0 X , i]/ y ) = \O x ,V^p y ) of two-mode 
bosonic states (here subscripts indicate particular vibra- 
tional modes; in what follows the subscripts will be omit- 
ted for given ordering of modes). To be more explicit, the 
realization of the transformation V can be expressed as 
the mapping in the extended Hilbert space of two bosonic 
modes and internal electronic levels in the following form: 



|0, V>') <E> \b) = NP\ b) B A \ip, 0) ® | a) 



(4) 



where the operators A and B represent a sequence of four 
types of unitary operations. The final projection fk\ on 
the state |6) selects conditionally the right outcome (A/" 
is a proper normalization constant). In the subsequent 
step one could use the two-mode linear coupler based 
on laser stimulated Raman transitions urn to swap the 
states of the vibrational modes, i.e., \0,ip^~ -> \ip',0). The 
additional 7r-pulse can be used to flip the electronic state 
from the level \b) into the initial level \a). 
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The operators A, B, appearing in the Eq.(||) in- 
dicate four essential steps which lead us subsequently to 
the desired transformations. 
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FIG. 1. The action of the operator A* (m) on the row with 
a fixed number m of phonons in the mode x. The population 
from the superposition ^ Vn : m\n) y is "swept" down into the 
vacuum \Q) y in the mode y. 

Step A. The operator A "spreads" the amplitudes c m 's 
of the component number states |n x ,0) ® \a) over the 
whole Hilbert space H x ®TL y ® Hi n so that the entangled 
state of the composed system becomes^] 



N,N 



V n , T . 



, Cm. G 



|m, n) (g) \a). 



(5) 



This task can be done by the method of 1-D quantum 
state synthesis proposed by Law and Eberly [ij. An im- 
portant tool represents also the photon-number depen- 
dent interaction H^ 1,m ^ considered by Kneer and Law ||] 
which enables to address the subspaces with a fixed num- 
ber m of phonons in the x direction. The operator A can 
be written as 



.4 



n e 



: 2cxp(i0 m )^( TO )^ r i/2 



(0) 



m=Q 



electronic level \a) to |&) by means of U^j™\ Then to 
"spread" \m 1 0) ® \b) we apply 1-D quantum state syn- 
thesis associated with the action of the operator A(m) to 

"prepare" the superposition J2n=o Vn,m\ m > n ) ® \b)- Af- 
ter that the electronic level |6) is flipped back to \a) via 



the action of U 



(l,m) 

-7r/2 exp(i</> m ) 



where 



[here (f) m 's represent proper 

phase factors which will be discussed later, see Eq.jTl])]. 

The appropriate interaction parameters for our 1-D 
quantum state synthesis can be found when we solve 
the inverse task which is given by the inverse transfor- 
mation: \m,0) (g> \b) = A^(m)^2 n V n , m \m,n) ® \b). The 
inverse task is based on " sweeping" down the probability 
from the component states of the given superposition in 
the y mode into the vacuum. Therefore the subscripts 
of the unitary operators (p — 2, 3) in (^) indicate 
that interaction parameters g p r = e l ^|<? p T| have to be 
chosen in such way that after the action of the tlfffi the 
amplitude (population) of the component state |$) be- 
comes equal to zero. The action of the operator A'(m) 
is shown schematically in Fig. [I]. We have applied the 
procedure introduced by Law and Eberly for synthesis 
of 1-D bosonic states in the straightforward way. There- 
fore we refer readers to the original paper 0] for other 
details (generalized Hamiltonians to operate beyond the 
Lamb-Dicke limit can be found in j[3|). Note that in 
our case we apply the state-synthesis procedure only in 
y direction and the resulting state (^) remains unknown. 

To resume, the action of the operator A encodes the 
matrix elements V n .m (multiplied with the amplitudes 
c m 's of unknown state) into rows of the 2D vibrational 
"lattice" of number states \m,n). In the next step an 
appropriate superposing of columns within the 2D vibra- 
tional number-state "lattice" is required. 

Step B. In the second step the operator B creates the 
state in which the amplitudes of the states |0, n) <S> \b) 
(n = 0, ...A) are proportional to c' n = J2 m =o^n,m c m . 
The state of the system \4j b ) = B\tp A ) after this synthesis 
step reads 



A(m) 



t>, (2) 



(3) 



\m,N;b)^ \m,N-l;c) 



.a 



(2) 



(•'!) 



\m,l;b) \m.Q;c) 



(7) 



The subscript of the unitary transformation Uglr in- 
dicates the required setting of the corresponding inter- 
action parameter g\T = \giT\e l ^ m . In other transforma- 
tions U$ (p — 2, 3) the subscripts denote steps in ID 
quantum-state synthesis as explained below. In other 
words, within a particular subspace with a fixed num- 
ber m of phonons in the x direction we flip from the 



Tn our synthesis algorithm we neglect off-resonant transi- 
tions between internal levels in H ( - 1,m \ Strictly speaking, the 
equality sign applies only in the limit x/|3i| — * 00 ■ 



N N 

i^ B } = 7^=10,^)^} + EE 



2 m „ \m,n) &> | a; 



The operator B can be written in the form 

rjv-i 



B = 



n ^( m ) 



u 



(1.JV) 

-17T/2 



where 



(8) 



(9) 



B(m) = fit 1 -" 1 *) 

£>{lll) Ly _ iar ctan(l/v / A'+l-m) L/ - l7, '/ 2ly -i7r/(2V^+T)' 

Here the subscripts of the unitary operators Ug^l indicate 
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again the proper choice of interaction parameters £f p T.fJ] 
The action of the operator Bona particular column 
with a given number n of phonons in the y mode is shown 
schematically in Fig. 0. Simultaneously, the operator B 
acts in the same way on "parallel" columns with differ- 
ent n. As illustrated in the upper part of Fig. [| the 
operator U_ff 2 ^ OT N = 5) transfers the population 
from the internal level |a) to \b) only in the row with the 
fixed number of quanta to = 5 in the mode x. Further, 
this population is transfered to the internal level \b) in 
the neighboring row with the number of quanta to = 4 
performing transformations U_) n , 2 U^ 4 ' 
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FIG. 2. The action of the operator B on the column with 
n vibrational quanta in the mode y. The same changes occur 
simultaneously on "parallel" columns with different n. The 
upper figure shows the action of U_^) 2 (for N — 5) followed 

by the basic sequence of operations B(m) (for m = 4). De- 
creasing the number of quanta m in x direction the basic 
sequence B (m) is recursively repeated as shown in the lower 
figure. After the action of the operator B the value of the 
amplitude s„ = Q2 m =o Vn )/VW+T. 



^Here and 



m 



Fig. 



the required setting g^r = 



— iir/(2\/m + 1) refers for clarity to the Lamb-Dicke regime 
rj x <S 1- Outside of the Lamb-Dicke regime the Rabi fre- 
quency \Jm + 1 is simply replaced by the nonlinear Rabi fre- 
quency e'^^im^J/v^nTT where L„ denotes the associ- 
ated Laguerre polynomial. 



The equal superposition of the amplitudes V n ^c^ and 
V n ^c^ in the row with m = 4 is obtained after ac- 
tion of i.e. undergoing one half of the 

-iarctan(l/\/2)' b to 

Rabi flipping. Decreasing m (the number of quanta in 
the mode x) the basic sequence B(m) is recursively re- 
peated as indicated in the lower part of Fig. @. Note 



that U 



(i,m) 



is responsible for adding of 



- i arctan( 1 / y/N-m+1) ' 

the amplitude V n ,m,c m (with a proper weight) to previ- 
ously superposed amplitudes Ylk=m+i Vn,kCk performing 
an adequate part of the Rabi flipping. After the ac- 
tion of the whole operator B the value of the amplitude 



-T N V 



c m . The transformed state \ip 



is given in Eq.(|8|). 

At this step we should notice that each action of the 
"elementary" unitary operator [A 1 '™), associated with 
the interaction Hamiltonian fl"( 1 > m ) with A y = —s m , 
causes significant phase shifts on off-resonant rows with 
the number of quanta in the mode x different from 
to ||. These phase shifts have to be compensated 
in advance in order to "superpose" the amplitudes via 
the action of B(m) as described above [see the role of 
U. m \ , nr, -=-,1. This compensation can be done 

— i arctan(l/ y/N— m+l) J r 

when we include appropriate phase shifts (f> m 's directly 
in the operator A [see Eq.(^)]. The explicit expression 
reads 



m-l N 



f 



k=0 



k—'m+l 



(11) 



where /mTlffil 7 ") = arg[cos(f2^ ; r 
,(*) 



2£2 



TTT 



sin(fl 



(fe). 



with Um = \J{sk — s m ) 2 /4 + | <7x | 2 - The origin of the ex- 
pression ( pi] ) can be traced back to the operators [A 1 '*"') 
in the steps A and B. The aim is to cancel the phase 
shifts of the amplitudes in order to "superpose" them on 
the m-th row by means of B(m). Therefore the first sum 
in (jll]) compensates (in advance) the subsequent shifts 
in A due to U^ 1 ' k > for k = 0, . . . , m — 1. The second sum 
compensates the shifts which will take place during "su- 
perposing" operations B(k) for k = m + 1, . . . , N which 
precede B(m). 

Step C. Comparing the state \ip B ) [Eq.(|])] with the de- 
sired one [Eq.(||)] we see that the target state is entan- 
gled to the internal level \b). On the other hand, also 
undesired component states |m 7^ 0, n) ® \a) are now 
contributing into Eq.(|8|) with nonvanishing amplitudes. 
However (fortunately), all the undesired components are 
entangled with the internal level \a). Therefore we can 
perform a conditional measurement to project the state 
vector (||) on the internal level \b). To be more specific, 
the internal state of the ion can be determined by driv- 
ing transition from the level \a) to an auxiliary level |r) 
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and observing the fluorescence signal^. No signal (no in- 
teraction with probing field) means that the undisturbed 
ion is occupying the level \b) being in the motional state 
lO, - ;//). It means that after the conditional measurement 
[indicated in (||) by the projector P^] the state vector (||) 
is reduced to the desired state vector (^) . The probability 
to find the right outcome for the unitary transformations 
is equal to j^rf- 

In spite of the involved conditional selection of the right 
outcomes, our algorithm is universal as the sequence of 
the "elementary" operations (with appropriate interac- 
tion parameters) which represents the desired transfor- 
mation is always independent of input states. Moreover, 
in contrary to conditional measurement schemes known 
from quantum state preparation, in our case the proba- 
bility of the right outcome is constant being also inde- 
pendent of input states. 

IV. DISCUSSION 

One of important applications of the operator synthe- 
sis is a realization of universal quantum gates for qubits 
which are encoded in vibrational levels. The number 
states of the vibrational mode can represent a quantum 
register. 

To illustrate our synthesis procedure we considered a 
realization of the operator which "rotates" the popula- 
tion between N vibrational Fock states under considera- 
tion. 

1 

V R = 7^=(E M + !>01 + |0><ATj). (12) 
ViV + 1 j=0 

It corresponds to a cyclic "rotation" of the quantum reg- 
ister. The operator Vr represents the unitary exponential 
phase operator of the Pegg-Barnett formalism . 

In the second example the synthesis procedure is ap- 
plied for the unitary op erator of the quantum Fourier 
transform defined as |T(| 

Vqft = TWti £ exp ( 27 %TI J |m)H (13) 

The operator of the quantum Fourier transform repre- 
sents an important tool in quantum computing jl6) . 




0.2 1 ■ ■ ■ 1 1 
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FIG. 3. The fidelity of outputs to the ideally transformed 
states as function of the ratio f° r the quantum Fourier 

transform Vqft (thick solid line) and the exponential phase 
operator Vr (thin solid line). Lamb-Dicke parameters are 
V* = Vy = 0-4. 

In the presented synthesis procedure we have neglected 
transitions between internal levels on off-resonant "rows" 
when the interaction Hamiltonian H^ 1 '" 1 ) is applied. 
Strictly speaking, the off-resonant transitions in H^ 1 *™ 1 ) 
can be neglected only in the limit — > 00. Our 
estimation of the error due to the finite values of the ra- 
tio x/lffi| (feasible in practice) is based on the fidelity of 
the outputs to the ideally transformed states (Q). The 
fidelity of two states |$), |$') is defined as their squared 
scalar product |(<i>'|$)| 2 . As a testing input state of the 
x-mode we can take a uniform superposition of involved 
number states, i.e., = J2m=o \ m )x (the initial 

internal state of the ion is \a) and the vibrational y-mode 
is in the vacuum). The figure ||| shows that the fidelity is 
close to one for large ratios x/lffil ^ 1 f° r both opera- 
tors Vr and Vqft|] In such case both the validity of the 
approximate Hamiltonian HO-'™) is justified S and the 
level flipping (|a) <-> \b)) on off-resonant "rows" (as the 
source of non-ideal fidelity) can be neglected. 

Let us note that for nonunitary transformations the 
probability depends on initial states and can be from the 
interval (0, 1). However, the realization of the transfor- 
mation (given by the sequence of the "elementary" oper- 
ations with appropriate interaction parameters) is always 
independent of input states. 



'As a check one could drive also transition from the level 
jc) to another auxiliary level \r'). Errors in the synthesis 
procedure are thus indicated by presence of the fluorescence 
signal. 



'In Fig. ^ we consider parameters outside of the Lamb-Dicke 
regime. To operate in the Lamb-Dicke regime requires a fur- 
ther increase of the ratio ^ 1 to reach a fidelity close 
to one. 
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V. CONCLUSIONS 

In this paper we have proposed a constructive algo- 
rithm for synthesis of operators, a superior task to syn- 
thesis of quantum states. Our method allows us to find 
analytical expressions for switching times and interaction 
parameters of the utilized laser stimulated processes via 
which an arbitrary unitary dynamics can be realized. 
One of important applications of the operator synthe- 
sis is a realization of universal quantum gates for qubits 
which are encoded in vibrational levels. As an example 
we consider realization of the discrete Fourier transform. 

The solution of the problem we present in our paper is 
neither unique nor optimal (in a sense of number of ele- 
mentary operations used for a construction of the given 
unitary operator) . The optimization of the procedure is 
the problem which has to be solved. The other problem 
which deserves attention is the stability of the algorithm 
with respect to noise inherent in the system. In fact, one 
can consider two types of uncertainties which might play 
an important role. Firstly, it is the noise induced by the 
environment, i.e. the elementary gates are not unitary. 
The second source of noise (kind of technical noise) is due 
to the fact that it is not possible to keep the interaction 
times and parameters to be fixed as given by the the- 
ory. Fluctuations in these parameters might reduce the 
fidelity of the realization of the desired unitary evolution. 
We will address these questions elsewhere. 
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